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It is shown that three-dimensional nonparaxial beams are described by the 
oblate spheroidal exact solutions of the Helmholtz equation. For the first time, 
their beam behaviour is investigated and their corresponding parameters are 
defined. Using the fact that the beam width of the family of paraxial Gaussian 
beams is described by an hyperbola, the connection between the physical 
parameters of nonparaxial spheroidal beam solutions and those of paraxial 
beams is formally stablished. These results are also helpful to investigate the 
exact vector nonparaxial beams. © 2008 Optical Society of America 
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Recent developments in science and technology brings the necessity of revisiting 
the theory and concepts of nonparaxiality of optical beams. Its understanding has 
great relevance in the description of optical fields which are tightly focused or, optical 
beams whose diameter can be of the order of a few wavelengths as can be the case in 
optical nanolasers. 

When studying beams beyond the paraxial approximation, the standard ap- 
proach, introduced by Lax, has been to start with the paraxial solution and then 
to include some corrections.^*^ Rigurously speaking, nonparaxial beams are solu- 
tions of the wave equation without the paraxial approximation, in other words, they 
are solutions of the Helmholtz wave equation V 2 E + k 2 E = 0, with k the wave num- 
ber. Thus, such kind of solutions must be investigated in order to describe nonparaxial 
beams. 

A class of nonparaxial solutions has been introduced assuming a point source 
located at a complex position along the z— axis. The spherical wave field from this 
virtual complex source is converted into a directional wave field which remains a 
rigurous solution of the Helmholtz equation.^ This virtual source mode has wave- 
fronts that are approximately oblate spheroidalPThis fact has also been discussed by 
several authors and even have expressed the virtual source mode wave in spheroidal 
coordinates) 7 * 8 1 9 1 ^ Whatever the coordinate system used, the field expression reduces, 
under the paraxial approximation, to Gaussian beam. ^ * 8 * 9 * 1 " * 11 * 

The virtual complex source point solutions carry an inherent singularity which 
makes them inadequate to describe propagating fields near the origin or focal source 
point To elliminate this problem a non-singular superposition of incoming and 

2 



outgoing spherical waves has been usedj^ however even this non-singular solution 
has problems since, to realise them physically, infinite energy is required.^ 

In this letter we investigate the solutions of the scalar Helmholtz equation in 
oblate spheroidal coordinates and we show that they can have a beam behaviour. 
By identifying that the evolution of the beam width of a Gaussian mode follows an 
hyperbola we relate the physical parameters of nonparaxial spheroidal beam solutions 
and those of paraxial beams. This also allows to define quantitatively a threshold 
between nonparaxial and paraxial beams. In the paraxial limit the oblate spheroidal 
solution tend to Laguerre- Gaussian beams. 

Pure plane, spherical or cylindrical waves of the Helmholz equation, cannot be 
used, in a strict sense, to describe optical beams if we understand as such those field 
distributions that are concentrated around an imaginary line that will be identified 
as the beam propagation axis. Here we show that solutions of the Helmholtz equation 
in spheroidal coordinates can satisfy this requirement depending on the value of the 
ellipsoidal parameter (to be defined below). 

The propagation axis will be the z— axis that will also be the azimuthal axis of the 
oblate spheroid coordinate system. This coordinates are defined by x = <i(^ 2 +l) 1 / 2 (l — 
r] 2 ) 1 / 2 cos0, y = d(£ 2 + 1) 1 / 2 (1 — 77 2 ) 1 ' 2 sin0 and z = d£r), where d > is the distance 
from the origin to the foci! i3 | 14 l The variation range of the coordinates is determined 
by the physical problem. In the present case they are < 77 < 1, — 00 < £ < 00 and 
< cf) < 2n. The first two coordinates can be related to an eccentric radial coordinate 
p = sinh -1 £ and an eccentric angular coordinate by 9 = cos -1 77. From the above 
coordinate transformations, the condition = 1 defines the z— axis. 
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In oblate spheroidal coordinates (£, 77, 0) the Helmholtz wave equation V 2 E + 
k 2 E = 0, is separated into the next set of equation d 13 * 14 * 

^| + m 2$ = o 

dtp 

where c = kd with a m „ is a constant parameter. The spheroidal parameter c is a 
meassure of how far is the spheroid from a sphere for which d = and in consequence 
c = 0. From Eq. (JTJ the solutions of Helmholtz equation have the form of a Lame 
productP' 

E nm (£, 77, 0) = fl« (£; c) S m „ (77; c) (2) 

where in the spheroidal radial functions R$ n (£; c), the superindex p defines the kind 
of solution (first kind p = 1, second kind p = 2 and so on). With p = 3, 4 the solutions 
are linear combinations of the first two, namely R^^ 1 = RmL ^ ^EL anc ^ t nese are 
what we will use as they, can be used to describe travelling waves. There are two 
spheroidal angular solutions (r); c) but, we exclude the second one since it tends 
to infinity as it approaches the z-axis then, only that of first kind can describe beams 
of finite energy and so we drop the superindex for this solution. 

In order to understand the physics of Eq. (J2J) we refer to Fig. 1 where the intensity 
of the modes with m = n = and m = n = 1 are shown. One immediately realises 
that for c 7^ 0, the wavefield distribution is concentrated around the z— axis, in other 
words, it behaves like a beam. The angular part of (J2J) determines the transverse 
pattern of the wavefield while the radial part defines the propagation features. The 
azimuthal and radial parts of the solution determine the form of the propagating 
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phase front, when m = they are oblate spheroidal, Fig la), and when m 7^ the 
solution describes wavefields with rotating wavefronts, Fig lb). 

To relate spheroidal solutions, Eq. (J2J) with those of the paraxial wave equation 
analyse the geometry of paraxial beams. Their beam width w(z) is given by w 2 (z) = 
Wq(1+z 2 /Rq) where wo is the beam waist and Rq is the Rayleigh distance.^ Rewriting 
this expression we get the hyperbola equation w 2 (z)/wq — z 2 /R\ = 1 to which we 
can associate an oblate spheroidal coordinate system, this is shown in Fig. 2. In it, 
d is the distance to one of the foci and w is the vertex of the hyperbola. From the 
geometry of the hyperbola, the equation 

R 2 = d 2 - w 2 (3) 

has to be fulfilled. From the equation of the hyperbola, the angle of the asymptote 
or far field angle is tana = wq/Rq, this assigns a geometric place to Rq in Fig. 2. 
We now make the connection between paraxial and spheroidal beams through the 
common parameter d. Sustituting c = kd and Rq = kw\/2 in Eq. 0, after some 
algebra, we get 

y/2(-l + VT+di) 
w = J . (4) 

This equation is the link that connects oblate spheroidal beams and paraxial beams. 
There must exist a range of values of c for which oblate spheroidal beams describe 
paraxial beams. This can be obtained, expressing this parameter in terms of the far 
field angle c = 2/(tanasino!), since paraxial beams cannot diverge at cone angles 
larger than a ~ 0.5rad ~ n/$^we have c p ~ 7 this value defines the frontier between 
paraxial and nonparaxial beams. In Fig. 3 we show the plots of the on-axis intensity 
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and transverse profile S mn (77; c p ) of an oblate spheroidal beam in this 
limit and compare them with the corresponding paraxial beam. Also, for reference, 
it is plotted the normalised axial intensity for a highly nonparaxial beam with c < c p 
(inner thin line). The small differences can be attributed to that the paraxial Laguerre- 
Gauss beam is in the limit of its validity. In fact, for c > c p the behaviour of the oblate 
spheroidal beams is even closer to paraxial beams. We have evaluated the wavefield 
profiles for large values of c > c p and found that the transverse part of the spheroidal 
beams described by the angular solutions S mn (77; c) tend to the Laguerre- Gauss beam, 
as expected.^ This similarity occurs earlier for the zero-order mode for c just above 
c p . Finally, notice that from Eq. (J3J) when c — > 00 the radial functions tend to 

plane waves for the fundamental order mode and when c — > 0, the nonparaxial oblate 
spheroidal beams tend to resemble spherical waves. 

That the solutions discussed here are different from those obtained with the 
complex source point method (CSPM) expressed in spheroidal coordinates.'^^^'The 
zero-order nonparaxial beams obtained with the CSPM have a singularity inherited 
from its construction from a source point The oblate spheroidal beams presented 
here do not present any singularity whenever c ^ 0. Using the CSPMjZEEEEl h j s 
usually proposed that in the paraxial limit d can take the value of Rq. From eq. (j2J) 
it is clear that this situation can never be possible since it would imply a physical 
inconsistency, w = and a far- field angle a = (see Fig. 2). 

In conlcusion, for the first time, we have demonstrated that spheroidal beams 
describe paraxial and nonparaxial beam exact solutions of the Helmholtz equation. 



i? (3 ' 4) f^-c ) 
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We have established a threshold that separates nonparaxial and paraxial regimes. 
Higher oblate spheroidal modes reduce in a natural way to Laguerre-Gauss beams 
when the spheroidal parameter c exceeds the critical value c p . The results presented 
here were in the neighborhood of the paraxial limit and thus the scalar treatment is 
valid. Our approach extends straightforward to vector nonparaxial beams in Eq. (j2J, 
to obtain the corresponding vectorial components obtained.^ Thorough investigations 
are underway. 

Authors acknowledge to CONACyT. GRM Acknowledge support to CONACyT 
grant 118631/120329. 
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Figure Captions 

1. Intesity and phase evolution (left), and transverse pattern (right) of oblate 
spheroidal modes for c = 7; a) m — n — b) m — n — 1. 

2. Hyperbolic-elliptic geometry of a paraxial beam defined by the beam width 
evolution. 

3. a) Comparison of on- axis intensity for oblate spheroidal nonparaxial (thick solid 
line) and paraxial (dashed line) beams. Also it is shown the normalised intensity 
for a highly nonparaxial beam (thin line), b) Corresponding transverse ampli- 
tude profiles at the plane z — 0. Notice that for the Laguerre paraxial mode r 
can be obtained from r 2 = x 2 + y 2 = d 2 (1 — if). 
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